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5t , We use the recently derived CP phase dependent analytic results for the super- 

symmetric electro- weak correction to — 2 to constrain the explicit CP phases in 
softly broken supersymmetry using the new physics effect seen in the g-2 Brookhaven 
measurement. It is shown that the BNL data strongly constrains the CP violating 
phase 9fj_ (the phase of the Higgs mixing parameter ^u) and ^2 (the phase of the 
SU(2) gaugino mass 1712) and as much as 60 — 90% of the region in the ^2 — ^/^ plane 
is eliminated over a significant region of the MSSM parameter space by the BNL 
constraint. The region of CP phases not excluded by the BNL experiment allows 
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for large phases and for a satisfaction of the EDM constraints via the cancellation 
mechanism. We find several models with large CP violation which satisfy the EDM 
constraint via the cancellation mechanism and produce an a^^^ consistent with 
the new physics signal seen by the Brookhaven experiment. The sparticle spectrum 
of these models lies within reach of the planned accelerator experiments. 

1 Introduction 

Recently the Brookhaven experiment E821 made a precise determination of the muon 
anomaly = {g^^ — 2)/2 and found a deviation from the Standard Model result at the 
2.6a level. The experiment finds|I[] 

A correction to is expected in supersymmetric models and a realistic analysis of the 
correction to was given in the supergravity unified model with gravity mediated break- 
ing of supersymmetry in Refs.0, Specifically in the analysis of Ref.0] it was pointed 
out that the supersymmetric electro-weak correction to could be as large or larger 
than the Standard Model electro-weak correction The fact that the supersymmetric 
electro-weak contribution to can be large is supported by several later works which 
include constraints of the unification of the gauge coupling constants using the high pre- 
cision LEP data, radiative breaking of the electro-weak symmetry, experimental bounds 
on sparticle masses and relic density constraints on neutralino dark matter assuming R 
parity conservation 0, §]. Soon after the release of data[|l| on the observation of a differ- 
ence between the experimental value and the theoretical prediction of in the Standard 
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Model several analyses appeared [9-17] both within the super symmetric framework [9-15] 
as well as in non-supersymmetric scenarios [|l^, [l^. Specifically the work of Ref. [^ used a 



2(T error corridor on a^"*^ where a^"*^ = a^^'P — a^^'^ to constrain SUSY, i.e., 



10.6 X 10"^° < aj^'y < 76.2 x 10"^° (2) 



where the error corridor also includes uncertainties due to hadronic error []TB| in the the- 
oretical predictions. The analysis of Ref. [^ showed the interesting result that the BNL 
data implies that the sparticles have upper limits which lie within reach of the planned 
accelerator experiments. Thus, for example, one finds that within the minimal SUGRA 
model]!^ the BNL constraint implies that the lighter chargino mass < 600 GeV, 



mi < 750 GeV and mo < 1.1 TeV (tan f3 < 30) consistent with the fine tuning criteria ||20 
These mass ranges are within reach of the LHC and a part of the parameter space may also 
be accessible at RUNII of the Tevatron. Further, it was shown that the BNL data implies 
that sign{fi) is positive using the standard sign convention [pi|]. This result is consistent 
with the experimental 6 — s + 7 constraint that eliminates much of the parameter space 



for the case when yU < 0||2^ . It was also discussed in Ref.|P that the effects from extra 



dimensions ||23|| on — 2 are typically very small and do not pose a serious background 
to the supersymmetric electro-weak contribution to — 2. 

The analysis carried out in Ref. [0 was under the assumption of CP conservation where 
the phases of all the soft SUSY parameters are set to zero. However, in general the soft 
breaking parameters can be complex and their presence brings in new sources of CP 
violation over and above the one in the CKM matrix of the Standard Model. The normal 
size of the phases is 0(1) which creates a problem in that phases of this size typically lead 



to the electric dipole moment for the electron and for the neutron which are in excess of 
the experimental limits[^]. Some of the possible ways to avert this disaster consist of 
assuming small phasesp5|, assuming that the sparticle spectrum is heavy p6[| and more 



recently the possibility that there are cancellations! 27] which allows for large phases and 
a not too heavy sparticle mass spectrum. For the last scenario one will have then large 
CP violating phases which would affect a variety of low energy physics, such as sparticle 
masses and decays, Higgs mixing, proton decay, B physics and baryogenesis. Thus the 
inclusion of CP phases is an important new ingredient in SUSY phenomenology. Now 
a^^^^ also depends on the CP phases and thus the experimental constraints on a^^^^ 
can be translated into constraints on the CP phases. 

A full analysis of the effects of CP phases on (7^ — 2 in the minimal N=l supergravity 
and in MSSM was given in Ref.|Q. Remarkably it is found that the CP phases strongly 
affect a^^^^ in that they can change both its sign and its magnitude. It is this fact, i.e., 
that the a^J^^^ is a very sensitive function of the CP phases that leads us to utilize the 
current BNL data to constrain the CP phases. The phases that enter most dominantly in 
the Qfj, — 2 analysis are and ^2 and one finds that the BNL constraint eliminates a big 
chunk of the parameter space in the ^2 — O/^i plane. The domains allowed and disallowed 
by the BNL data depend sensitively on mo, mi and tan/5 and less sensitively on other 
parameters (Here mo is the universal soft breaking mass for the scalar fields, mi is the 
universal gaugino mass, and tan/5 =< H2 > / < Hi > where H2 gives mass to the up 
quark and Hi gives mass to the down quark and the lepton). Often as much 60 — 90% of 
the area in the ^2 ~ plane is excluded by the BNL constraint. In the limit of vanishing 
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phases the allowed region reduces to the constraint jj, > which was deduced in the earlier 
analysis under the constraint of CP conservation. Of course, not all the parameter space 
allowed by the — 2 constraint is allowed by the EDM constraints. However, we find 
that these constraints can be mutually consistent. Thus we give examples of models where 
the phases are large, ie., 0(1), the EDM constrains are satisfied and the contribution to 
Q^susy jg consistent with the new physics signal seen in the precise BNL measurement. The 
outline of the rest of the paper is as follows: In Sec. 2 we give some of the basic formulae 
which enter into the g^ — 'i analysis with CP violation. In Sec. 3 we give a discussion of the 
numerical results obtained by imposing the constraint of Eq.(2). Conclusions are given in 
Sec.4 



2 — 2 with CP Violating Phases 

Q^susy Qj^g loop level arises from the chargino exchange and from the neutralino ex- 

change. The chargino contribution is typically the dominant one although the neutralino 
contribution can become very significant in certain regions of the parameter space. In 
fact the neutralino exchange contribution is central in determining the boundary of the 
allowed and the disallowed region in the plane of the CP violating phases on which a^"*^ 
sensitively depends. To define notation and explain the main features of the analysis we 
exhibit below the CP phase dependent chargino contribution and the reader is referred to 



the Ref. |^ for the full analytic analysis including the neutralino exchange contribution. 



The chargino mass matrix with CP phases is given by 
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I m2 1 e*^2 v^miy sin (3 \ 
Mc=\ (3) 
v2mvi/cos/3 |/i|e*^f / 

where we have absorbed the phases of the Higgs sector by field redefinitions. In Eq.(3) 
6^ is the phase of the Higgs mixing parameter /x, ^2 is the phase of the SU(2) gaugino 
mass 7712- The chargino mass matrix can be diagonahzed by the biunitary transforma- 
tion U*McV^^ = diag{rh^+,rh^+) where U and V are unitary matrices. The chargino 
contribution is given by||28|| 

= a?" + (4) 



where 



a 



+ 



7W i=l 



and 



a 



n*- ^J_n, 1=-"-*^' 



Here F3(x) = (x - l)~^(3x^-4x+l-2x2/nx), F4(x) = (x - 1)"^ (2x^+3x^-6x+l-6x2/nx) 
and = m^/A/2Mvi/ cos /?. As discussed in Ref. ||28[| the entire phase dependence of the 



chargino contribution to resides in the combination 9^+^2- The neutrahno contribution, 
however, depends on additional phases and one can choose these additional combinations 
to be + ^1 and 6^ + a^o where ^1 is the phase of the U(l) gaugino mass mi and 
is the phase of the universal parameter Aq of trilinear soft SUSY breaking term in the 
scalar potential. 
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3 Discussion of Results 

We want to analyze the effect of the constraint of Eq.(2) on the CP phases on which 
a^^^^ sensitively depends. We shall work in the region of the parameter space where 
sparticle masses are of moderate size and CP phases are 0(1). In this region one can 
manufacture an a^^^^ of the size of the new physics effect indicated by the BNL experi- 
ment. However, since the sparticle masses are moderate size and the phases are 0(1) we 
need the cancellation mechanism to achieve consistency with the EDM constraints. For 
the purpose of the analysis we shall use the parameter space discussed in Ref. p9| which 
consists of the parameters: mg, ^1/2, Aq, tan/?, 6^, ^1, ^2 and ^3 where ^1 is phase of the 
U(l) gaugino mass mi, and ^3 is the phase of SU(3) gaugino mass fh^. The electro-weak 
sector of the model does not involve the SUSY QCD phase ^3 which, however, enters in 
the EDM analysis of the neutron. The neutralino exchange contribution depends also 
on ^1 and a^o in addition to its dependence on 6*^ and ^2- However, the dependence of 
the sum of both contributions on ^1 and is weak. Thus mainly the phases strongly 
constrained by the Brookhaven experiment are ^2 and 9^. 

In Fig.l we display the allowed parameter space in the ^2 — plane in the range 
— < ^2 ^ and — TT < 6^ < n for the specific input values of mo, mi/2, tan/5, Aq, 
and .^1 as given in the caption of Fig.l. As discussed in Sec. 2, the chargino exchange 
contribution to a^^^^ is a function only of the combination 9 ^+^2- This means that in the 
part of the parameter space where the chargino contribution is dominant a value of ^^ + ^^2 
allowed by the BNL constraint will generate a 135" line in the .^2 — 9^j_ plane. Similarly a 
range of allowed values of 9^ + ^2 will generate an area at 135° incline and we see that 
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is approximately true in Fig.l. Now, of course, if the chargino contribution was the sole 
contribution in a^^^^ Fig.l would consist of only parallel lines at 135*^ incline within 
the allowed range of + ^2 consistent with Eq.(2). However, a^^"^^ also contains the 
neutralino exchange contribution which is strongly dependent on ^2 and individually 
even when 0^ + ^2 is fixed. The strong dependence of the neutralino contribution on ,^2 
when 9^ + ,^2 is fixed is shown in Fig. 2. Because of this the sum of the chargino and 
the neutralino contributions does not possess the simple dependence on 6^ and ^2 in the 
sum form. Thus in Fig.l the boundaries at 135° are not exactly straight lines since near 
the boundary the neutralino contribution can move a^^^"^ in or out of the allowed range 
admitting or eliminating that point in the parameter space of the admissible set. Further, 
since the neutrahno contribution violates the simple dependence on 6^ + ^2 it destroys 
the translational invariance of a^J^^^ on 9^ (with 6^ + ^2 fixed). We see this violation in 
Fig.l from the asymmetrical endings of the allowed region, i.e., the lower right hand and 
the upper left hand of the admissible region, are not mirror reflections of each other. In 
addition to 6^ and ,^25 the parameters tan/3, mi/2 and mo also have a strong effect on a^. 
We study the effect of changes in these below. 

To study the effect of the dependence on tan (3 we carry out an analysis in Fig. 3 similar 
to that of Fig.(l) but with tan/? = 5 and with all other parameters fixed at their values 
in Fig.l. Now as pointed out in Refs.|^, |] a^^^^ has a strong dependence on tan/?. As 
shown in the first paper of Ref. this dependence arises from the chiral interference term 
in the chargino exchange contribution which is proportional to tan j3 for large tan (3. Thus 
a reduction in the value of tan j3 reduces the magnitude of af^"^^ and its relative smallness 
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results in a smaller range in 6'^ + ^2 around 6'^ + ^2 = consistent with Eq.(2). Further, 
since the magnitude of a^^^^ is smaller and closer to the lower limit of Eq.(2) it is more 
sensitive to the neutralino exchange contributions which can move it out of the allowed 
region more easily reducing the allowed region in the ^2 — pl^ine which is what we see 
in Fig. 3. 

Next we study the effect of changing the value of mi/2- An increase in the value of 
mi/2 increases the chargino mass and the neutralino mass which reduces the magnitude of 
a^^^^ . The reduction in the magnitude of a^^^^ leads to a smaller range for 6'^ + ^2 which 
is what we see in Fig. 4 relative to Fig.l. Finally we look at the effect of the variation of 
mo on a^^^^ ■ Now similar to the effect of increasing mi/2, an increase in the value of 
mo decreases the overall magnitude of a^^^^ bringing it closer to the BNL lower limit as 
given by Eq.(2). As in the analysis of Figs. 3 and Fig.4 the fact that the overall magnitude 
of a^^^^ is smaller means that changes in the phase can more easily move its value out 
of the BNL admissible domain thus reducing the allowed range of 6'^ + ^2- This is what 
we see in Fig. 5 where mo = 400 GeV. The same argument would indicate that a further 
increase in the value of mg should further decrease the allowed range of 9^ + ^2 which 
is what we find in Fig. 6 where mo = 600 GeV. We note that a^J^^^ is more sensitive to 
changes in mi than in tjiq. This was seen already in the analysis of Ref. where the 
upper limit of mo was found to be significantly larger than the upper limit on mi. The 
above explains why the allowed area in Fig.4 is smaller than in Figs. 5 and 6. This is so 
because in Fig.4 mi is close to its upper limit while is Figs. 4 and 5 mo is significantly 
lower than its upper limit. We notice also that the allowed regions in Figs. 5 and 6 
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consist of complete straight lines which means that the neutralino contribution role here 
is suppressed severely by increasing mo. 

Now not all the parameter space admissible by the — 2 constraint in Figs. 1-6 is 
admissible by the constraints on the electron and on the neutron EDM. To satisfy the 
EDM constraints via the cancellation mechanism we have to utilize also the parameter ^3 
along with 9^, ^2 and other soft parameters. (We include the two loop effects of Ref. in 
the analysis.). In Table 1 we exhibit five points (a-e) that lie in each of the allowed regions 
of Fig.l and Figs. 3-6, i.e., point (a) lies in the allowed region of Fig. 1, point (b) lies in 
the allowed region of Fig. 3 etc, which satisfy the EDM constraints and the corresponding 
value of a^^^^ lies in the BNL range of Eq.(2). Thus Table 1 gives five models which 
have large CP violating phases, their contributions to the EDM of the electron and of 
the neutron lie within experimental limits using the cancellation mechanism and they 
produce a SUSY contribution to — 2 consistent with the signal observed by the BNL 
experiment. The sparticle spectrum corresponding to cases (a-e) of Table 1 is shown in 
Table 2. One finds that in all the cases the sparticle spectrum is low enough that some if 
not all of the sparticles must become visible at the LHC, and in some cases the sparticle 
spectrum is low enough to even lie within reach of RUNII of the Tevatron. We note that 
using points of Table 1 one can generate trajectories using the scaling technique given in 



Ref. IP If where the EDM constraints are satisfied and thus one can use this technique to 
produce many more models of the type discussed above. 
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Table 1: Cases where the EDM and the g-2 experiments are satisfied 



(case) ^2, dfi, ^3 (radian) 


4, dn (ecm) 


„SUSY 


(a) -.63, .3, .37 


-4.2 X 10-27^ -5.3 X 10"26 


47.0 X 10"^° 


(b)-.85 ,.4 ,.37 


4.2 X 10^27^ 4 g X 10-26 


10.8 X 10-10 


(c)-.8 ,.2 ,1.3 


4.0 X 10-27^ 5.4 X 10-26 


12.2 X 10-10 


(d)-.32 ,.3 ,-.28 


-1.2 X 10-27^ 3.3 X 10-26 


20.1 X 10-10 


(e)-.5 ,.49 ,-.5 


1.8 X 10-27^ -6.6 X 10-27 


12.7 X 10-10 



4 Conclusion 

In this paper we have used the new physics signal seen by the Brookhaven g-2 experiment 
and the recently derived CP phase dependent analytic results on the supersymmetric 
electro-weak correction to (7^ — 2 to put limts on the explicit CP violating phases that 
arise from the soft SUSY breaking sector of MSSM. Using a 2 a error corridor around 
the observed effect we find that the BNL constraint excludes a large region in the ^2 — dfi 
plane. The amount of the region excluded depends sensitively of tan/?, mo, and mi and 
less sensitively on the remaining parameters. In most of the parameter space the excluded 
region is as much as 60 — 90% of the total area, i.e., the area mapped by — vr < .^2 < tt, 
—71 < Ofj, < n. In the limit when the phases vanish the allowed region limits to the case 
fi > deduced in the previous analysis using CP conservation p. We also show that the 
regions allowed by the BNL constraint contains points where the cancellation mechanism 
operates and provide examples of models with large CP violation consistent with the ex- 
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Table 2: Sparticle masses (in GeV) for cases (a-e) in Table 1. 



(case) 


Xl) X21 X31 X4 


Xl, xt 


P'l,P'2 


Ml, U2 


(a) 


98.2, 186.9, 389.8, 403.5 


190.2, 405.6 


145.0, 209.6 


628.5, 647.4 


(b) 


97.2, 184.0, 408.3, 426.6 


187.0, 426.6 


144.6, 209.1 


628.6, 647.5 


(c) 


213.8,421.0, 845.8, 852.2 


429.3,853.2 


232.0, 397.2 


1335.3, 1378.0 


(d) 


98.1, 186.0, 378.4, 393.4 


189.2, 395.3 


413.5, 440.3 


738.3, 754.4 


(e) 


98.3,187.1,393.7, 407.3 


190.4, 409.3 


609.1, 627.6 


863.2,877.0 



perimental EDM limits of the electron and the neutron and with SUSY contributions to 
consistent with the new physics signal seen by the BNL experiment. These models 
also possess sparticle spectra which lie within reach of collider experiments planned for 
the near future. 

This research was supported in part by the NSF grant PHY-9901057. 
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Figure 1: A plot of the allowed region (shaded) in the ^2 — plane allowed by the 
constraint of Eq.(2) when mo = 100, mi = 246, tan/5 = 20, Aq=1, = .3, and = -5 
where all masses are in GeV. 
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Figure 2: A plot of the chargino contribution (dotted line), neutralino contribution 
(dashed line) and the total a^^^^ (solid line) as a function of ^2 in the range — tt < ^2 < tt 
when 0^-\- ^2 — —1, friQ — 100, mi = 246, tan/5 = 20, = 1, ^1 = .4, a^o — -5, where 
all masses are in GeV. The small fluctuation of the chargino contribution from exact 
constancy is due to small rounding off errors in the numerical integration program. 
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Figure 3: A plot of the allowed region (shaded) in the ^2 — 0^ plane allowed by the 
constraint of Eq.(2) with all the same parameters as in Fig.l except that tan/3 = 5. 
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Figure 4: A plot of the allowed region (shaded) in the ^2 — ^fj, plane allowed by the 
constraint of Eq.(2) with all the same parameters as in Fig.l except that mi = 527 GeV. 
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Figure 5: A plot of the allowed region (shaded) in the ^2 — plane allowed by the 
constraint of Eq.(2) with all the same parameters as in Fig.l except that mo = 400 GeV. 
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gure 6: A plot of the allowed region (shaded) in the ^2 — plane allowed by the 
constraint of Eq.(2) with all the same parameters as in Fig.l except that mo = 600 GeV. 



22 



